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Neutrino-neutrino interactions can lead to collective flavour conversion effects in supernovae and 
in the early universe. We demonstrate that the case of "bipolar" oscillations, where a dense gas 
of neutrinos and antineutrinos in equal numbers completely converts from one flavour to another 
even if the mixing angle is small, is equivalent to a pendulum in flavour space. Bipolar flavour 
conversion corresponds to the swinging of the pendulum, which begins in an unstable upright position 
(the initial flavour), and passes through momentarily the vertically downward position (the other 
flavour) in the course of its motion. The time scale to complete one cycle of oscillation depends 
logarithmically on the vacuum mixing angle. Likewise, the presence of an ordinary medium can 
be shown analytically to contribute to a logarithmic increase in the bipolar conversion period. 
We further find that a more complex (and realistic) system of unequal numbers of neutrinos and 
antineutrinos is analogous to a spinning top subject to a torque. This analogy easily explains how 
such a system can oscillate in both the bipolar and the synchronised mode, depending on the neutrino 
density and the size of the neutrino-antineutrino asymmetry. Our simple model applies strictly only 
to isotropic neutrino gasses. In more general cases, and especially for neutrinos streaming from a 
supernova core, different modes couple to each other with unequal strength, an effect that can lead to 
kinematical decoherence in flavour space rather than collective oscillations. The exact circumstances 
under which collective oscillations occur in non-isotropic media remain to be understood. 

PACS numbers: 14.60.Pq, 97.60.Bw 



I. INTRODUCTION 

The "refractive index" caused by a medium has a 
strong impact on neutrino flavour oscillations 1, 2, 3. 
0, H, Q. This matter effect is a standard ingredient 
for neutrino oscillations in laboratory experiments and 
in astrophysical environments. However, when neutri- 
nos themselves form a significant "background medium" 
as in the early universe or in core-collapse supernovae, 
the oscillation equations become nonlinear, sometimes 
resulting in surprising collective phenomena. Based on 
Pantaleone's key observation that neutrinos as a back- 
ground medium produce a flavour off-diagonal refractive 
index Q , the behaviour of dense neutrino gases was in- 
vestigated in a series of papers by Samuel, Kostelecky 
and Pantaleone @, |, M, M, E3, M, EI M, GJ]- Tw0 
classes of collective effects were identified in these pa- 
pers: we shall call them "synchronised" and "bipolar" 
oscillations, respectively. 

Synchronised oscillations occur when the neutrino- 
neutrino interaction potential is large compared to the or- 
dinary oscillation frequencies in vacuum or in a medium, 
and a sufficiently large asymmetry exists between the 
neutrino and antineutrino distributions. As a result, 
all neutrinos and antineutrinos oscillate with the same 
frequency that is a certain average of the ordinary os- 
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cillation frequencies. In the spin-precession analogy of 
flavour oscillations, the flavour polarisation vectors of all 
neutrino modes form one big spin that precesses in a 
weak "external magnetic field." This big spin is held 
together by the strong "internal magnetic field" formed 
in flavour space by the strong neutrino-neutrino interac- 
tion [T3|. Various qualitative and quantitative aspects 
of synchronised oscillations and applications to neutrino 
flavour oscillations in the early universe [H, [lj| H(J, Hl[ 
were studied a few years ago, motivated by the question 
if neutrinos with chemical potentials reach flavour equi- 
librium in the early universe before the epoch of big-bang 
nucleosynthesis. 

Much less attention has been paid to bipolar flavour 
conversions, but very recently it has been recognised that 
this peculiar phenomenon likely plays a crucial role for 
supernova neutrino oscillations 0, 0, [24[ • In the sim- 
plest case, bipolar oscillations occur in a dense gas of 
equal numbers of neutrinos and antineutrinos of the same 
flavour. For a suitable mass hierarchy, even a small mix- 
ing angle can cause a complete conversion of both neutri- 
nos and antineutrinos to the other flavour. We stress that 
a non-vanishing vacuum mixing angle is pivotal; bipolar 
oscillations, although a quasi self-induced effect, will not 
occur if the vacuum mixing angle strictly vanishes. 

The relevant conditions for bipolar conversions proba- 
bly occur for neutrinos streaming off a collapsed super- 
nova core where one expects a hierarchy of average ener- 
gies (E Ve ) < (Ep B ) < (E Vx ) with V x = V^V^Vr^Vr [H]. 

Assuming equal luminosities for all species, the number 
flux of, say, and is each smaller than that of v e 
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and v e respectively. For the inverted mass hierarchy one 
would then expect bipolar oscillations driven by the "at- 
mospheric" neutrino mass difference and the small mix- 
ing angle O^. 1 

Bipolar flavour conversion occurs when the neutrino- 
neutrino interaction energy, fi = ^f^Gyriy, exceeds a typ- 
ical vacuum oscillation frequency u> = Am 2 /2E. Fur- 
thermore, the neutrino-antineutrino asymmetry must not 
be too large in a sense to be quantified later. Once these 
conditions are met, bipolar oscillations take place for 
astonishingly small values of the mixing angle and are 
nearly unaffected by the presence of an ordinary back- 
ground medium, even if it is much denser than the neu- 
trino gas [22j. These oscillations are also unrelated to a 
Mikheyev-Smirnov-Wolfenstein (MSW) resonance of the 
ordinary background medium. 2 

The notion of "bipolar oscillations" originates from the 
claimed numerical observation that the flavour polarisa- 
tion vectors of all neutrinos and of all antineutrinos align 
together to form two "block spins" that evolve separately. 
Indeed, the analytic descriptions of Refs. (l2l . Il4j were 
based on the study of a system of one neutrino and one 
antineutrino polarisation vector that were taken to rep- 
resent, respectively, the complete ensemble of neutrinos 
and antineutrinos. We find that this description of bipo- 
lar oscillations is incorrect and actually has never been 
explicitly demonstrated in the literature. In fact, each 
mode of the neutrino and antineutrino ensemble evolves 
differently and the term "bipolar" is a misnomer. On the 
other hand, the behaviour is bipolar in the sense that 
neutrinos and antineutrinos oscillate in "opposite direc- 
tions" and thus form two separate cohorts, even if they 
do not form two block spins. Therefore, we use the term 
"bipolar" to describe this collective phenomenon. 

The newly recognised role of bipolar oscillations rep- 
resents a change of paradigm for supernova neutrino os- 
cillations. Inspired by these exciting developments we 
turn to an analytic study of bipolar oscillations with the 



1 We concentrate here on 13-mixing because and v T streaming 
off a supernova core have equal spectra, behave equally in ordi- 
nary matter, and are known to be nearly maximally mixed so 
that the 23-mixing is not relevant in this context. The "solar" 
mass difference is much smaller and the corresponding hierarchy 
is normal so that 12-mixing is expected to lead to less prominent 
effects at larger distances from the supernova core. In general, 
however, one would need to perform a three-flavour treatment. 

2 Of course, it has long been recognised that the high density of 
neutrinos near a supernova neutrino sphere could lead to refrac- 
tive effects comparable to those of the ordinary medium, and 
nonlinear effects could play an important role 26, 27, 28]. The 
conclusion of these early works appears to be that the neutrino- 
neutrino term causes a small shift of the oscillation parameters 
where an MSW resonance occurs and hence a small correction 
to the ordinary matter effect. Subsequent numerical studies of 
nonlinear effects in the supernova hot bubble region actually re- 
vealed flavour conversion for surprisingly small neutrino mass 
differences, but a connection to the bipolar oscillation mode was 
not made [29l, l3fjl] . 



aim of providing a simple qualitative and quantitative 
understanding of the salient features of this puzzling ef- 
fect. Analytic solutions for certain cases have already 
been provided in the literature [HI, HH. However, a judi- 
cious choice of variables enables us to write the equations 
of motion in the form of an ordinary pendulum. This pic- 
ture allows one to grasp the salient features of the bipolar 
phenomenon at a single glance. Moreover, it allows one 
to calculate explicitly the dependence of the bipolar os- 
cillation period on the vacuum mixing angle and on the 
density of an ordinary background medium. 

When the fluxes of neutrinos and antineutrinos are not 
equal (as expected for supernova neutrinos), the equa- 
tions for the flavour pendulum remain the same, but the 
asymmetric initial conditions imply the presence of an in- 
ner angular momentum (i.e., spin) of the pendulum: the 
system is equivalent to a spinning top subject to a torque. 
If the top is not spinning, we simply recover the motion of 
an ordinary spherical pendulum (i.e., bipolar behaviour 
of a symmetric v-v system). Otherwise the motion is 
more complicated. If the spin is sufficiently large, the 
top precesses in the force field exerting the torque (i.e., 
synchronised oscillations). If the spin is too small, the 
top wobbles or even completely turns over (i.e., bipolar 
behaviour of an asymmetric system). 

In the most general case, however, different neutrino 
modes have different energies and, in a non-isotropic 
medium such as that encountered by neutrinos streaming 
off a supernova core, couple to each other with different 
strengths ( "multi- angle case"). In this situation kine- 
matical decoherence rather than collective oscillations 
can obtain and is an unavoidable outcome certainly in 
the simplest system of equal numbers of neutrinos and 
antineutrinos. However, in the more realistic case of 
unequal neutrino and antineutrino fluxes in conjunction 
with a slowly varying effective neutrino density, collec- 
tive oscillations often still obtain. The exact criteria that 
determine if kinematical decoherence or collective oscil- 
lations occur remain to be understood. 

We begin in Sec. |TT] with the simplest bipolar system 
consisting of one polarisation vector for the neutrinos and 
one for the antineutrinos and establish the equivalence to 
a flavour pendulum. We then study analytically the im- 
pact of an ordinary background medium in Sec. IIIII In 
Sec. IIVI we allow the neutrino density to vary and show 
that this effect is crucial for the nearly complete flavour 
conversion in supernovae. In Sec. [V] we consider a system 
with an initial v-v asymmetry and show its equivalence 
to a spinning top subject to a torque. We apply our 
insights to explain the salient features of flavour conver- 
sion of the neutrinos streaming off a supernova core. In 
Sec. I VII we consider a system of many modes and discuss 
the conditions under which it is equivalent to the simple 
bipolar system. In Sec. IVIII we discuss the possibility of 
flavour conversion initiated by quantum fluctuations, in 
the absence of flavour mixing. A summary of our findings 
is given in Sec. IVIII1 
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II. BASIC BIPOLAR SYSTEM 
A. Equations of motion 

We begin with the simplest bipolar system initially 
composed of equal densities of pure v e and P e . All of 
them are taken to have equal energies so that the vac- 
uum oscillation frequencies are the same for all modes. 
We describe the flavour content of these ensembles with 
polarisation vectors in flavour space P and P, where over- 
barred quantities refer to antiparticles here and hence- 
forth. Without loss of generality we take these vectors 
to have unit length. As usual, the ^-component of the 
polarisation vector represents the flavour content of the 
ensemble, i.e., the survival probability of v e at time t is 
i [f + P z (t)]. We take the positive z direction to represent 
the electron flavour so that both P and P are initially 
unit vectors in the z-direction. 

In the absence of ordinary matter the general equations 
of motion Eq. (|A8j) are 

d t P = [+ojB + (i (P - P)] x P , 

8 t P = [-wB + a*(P-P)] xP, (1) 

where u > is the vacuum oscillation frequency, /i = 
\/2GFn u represents the strength of the v-v interac- 
tion, i.e., the density of the neutrino gas, and B = 
(sin 2^o, 0, — cos2#o) with vacuum mixing angle 9q. 

A mixing angle close to zero corresponds to the normal 
mass hierarchy in which v e is essentially identical with the 
lower mass state, while 9q near 7r/2 corresponds to the 
inverted hierarchy with v e residing largely in the heavier 
mass state. In the latter case we will also use the notation 

9 Q = tt/2 - O . (2) 

Therefore, whenever 9q appears as a small quantity it 
signifies directly that we are in an inverted mass situa- 
tion. Alternatively, one could restrict the vacuum mixing 
angle to < 6q < 7r/4 and switch to the inverted hier- 
archy with the replacement to — > —lu. However, in our 
treatment it is more natural to keep to always positive 
and extend the range of mixing angles to < #o < 7r/2. 

In order to illustrate the phenomena we wish to study 
we show in Fig. [1] the evolution of P z for a vacuum mix- 
ing angle near 7r/2, corresponding to an inverted mass 
hierarchy. At first P hardly moves at all, but after some 
time it flips almost completely. Therefore, even a very 
small mixing angle leads to complete flavour conversion. 
Of course, this simple system is periodic so that the mo- 
tion then reverses itself. The behaviour of P z is identical 
to that of P z . In other words, both v e and v e convert 
simultaneously to and v^. 

This evolution can be understood from Eq. ([TJ) . Ini- 
tially the difference of the polarisation vectors 

D = P P (3) 
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FIG. 1: Evolution of P z and P z for the system of equations 
Eq. fT} with #o = 0.01 (i.e., inverted hierarchy), u = 1, and 
strong neutrino-neutrino interaction /i = 10. (This figure is 
essentially identical to Fig. 3 of Ref. [171].') 



vanishes so that there is no neutrino- neutrino effect. 
When the polarisation vectors P and P begin to pre- 
cess in opposite directions around B, a non-zero D de- 
velops in the y-direction orthogonal to B. Both P and P 
then tilt around D, leading to a complete flavour rever- 
sal (inverted hierarchy) or to small oscillations (normal 
hierarchy) . 

Observe that P and P behave symmetrically and their 
z-components develop identically. This suggests that in- 
stead of these polarisation vectors one should use their 
sum and difference vectors as independent variables, i.e., 
D as defined in Eq. ©, and 3 

S = P + P . (4) 

The z-component of S quantifies the flavour content of 
the combined v e and v e ensemble. The equations of mo- 
tion for the S and D vectors are 

S = wBxD + fiDxS, 
D = uBxS. (5) 

The first line of Eq. ([5]) suggests yet another vector to 
describe the ensemble, 

Q = S -B. (6) 

M 



3 Our vector S corresponds to S— of Ref. [23| . whereas D corre- 
sponds to their S_|_, i.e., what Ref. 1221 calls a sum is what we call 
a difference, and vice versa. This reversal of roles arises because 
we work with polarisation vectors, while Ref. [22ll uses "neutrino 
flavour isospin (NFIS)" vectors. We discuss the advantages and 
disadvantages of these languages in Sec. I VIII 



4 



For strong neutrino- neutrino interactions (/i/w > 1) we 
can think of Q as identical to S. 

Since B = 0, it follows that S = Q. With B x Q = 
B x S, the equations of motion are now 



Q = fiT> x Q , 
D = uBxQ. 



(7) 



Clearly, the length of Q is conserved and stays at its 
initial value 



Q = |Q| = 



A* 



-I 1/2 



4 — cos 29 



(8) 



where we have used |B| 2 = 1, and the initial values |S| 2 = 
4 and B • S = -2cos26> n . 



B. Spherical pendulum 

The vector Q in flavour space plays the role of a spheri- 
cal pendulum in that its length is conserved so that it can 
move only on a sphere of radius Q. In this picture, the 
role of the different quantities is most easily understood 
if we consider the total energy of the system, 



H = wB-Q + -D 2 . 

2 



(9) 



up to a constant. The first term is the potential energy of 
the pendulum in a homogeneous force field represented 
by wB. The second term is the kinetic energy, with D 
playing the role of the pendulum's orbital angular mo- 
mentum. Observing that D • Q = — (u>//i)D • B is con- 
stant due to Eq. ([7]) and thus zero in our case D(0) = 0, 
the first line in Eq. ([7]) implies 



1 Q x Q 

(x Q 2 



(10) 



and hence D 2 = /i~ 2 Q 2 /Q 2 . The scale of the poten- 
tial energy is set by the vacuum neutrino oscillation fre- 
quency u>, whereas / = is to be identified with the 
moment of inertia. The latter should be compared with 
/ = mi 2 for an ordinary mass suspended by a string of 
length £. The role of inertia in the pendulum analogy is 
played by the inverse strength of neutrino-neutrino in- 
teraction! 



C. Plane pendulum 



Equation (jTTj) can be further simplified to 

i/3 = — k 2 sin((/> + 26*o) , 

where 

K 2 = LdflQ . 



(12) 



(13) 



The inverse of k is the characteristic time scale for 
the bipolar evolution. In the limit of strong neutrino- 
neutrino coupling, Q « 2 and hence k ss \J2lo[l. In the 
opposite limit (/i <C uj), we have Q ~ uj//i so that n sa u). 
In this latter case the characteristic frequency of the sys- 
tem is the vacuum oscillation frequency to, and P and P 
oscillate independently. 

The equations of motion ([IT]) follow directly from the 
classical Hamiltonian for a simple pendulum 



H{<p,D) 



-[1 



cos((/) + 20 o )] 



1 9 

2^ 



(14) 



where tp is a coordinate and D its canonically conju- 
gate momentum. Indeed, the equations pip are but 
<p = dH/dD and D = -dH/dip. 

Assuming a small vacuum mixing angle 9q and a small 
excursion angle ip of the pendulum, the potential can be 
expanded 

,2 



V(<p) 



ir [1 

„2 



cos((/> + 29 )\ 
+ 29 Q f + ... ■ 



(15) 



In this case the system is equivalent to a harmonic os- 
cillator with frequency k. On the other hand, for angles 
near it so that |</> + 28q — 7r| <^ 1, we get the same expan- 
sion but with a negative sign; the system corresponds to 
an inverted harmonic oscillator. 



D. Bipolar flavour conversion 

1. Normal hierarchy 

Consider the small mixing limit in the normal hierar- 
chy. Here, the initial condition </>(0) « —(uj/^iQ) 29q > 
—29q puts ip near the minimum of the potential V(ip) at 
t = 0. Since ip(0) = 0, the system remains trapped inside 
the cosine potential well, oscillating about the minimum 
fmin = —29 with amplitude (1 — uj/(j,Q) 29q and fre- 
quency k. In terms of P z and P z , any departure from the 
initial P z = P z = 1 is at most second order in 9$. No dip 
features develop in this scenario. 



Our initial conditions P(0) = P(0) = (0,0,1) imply 
D(0) = 0. The pendulum's subsequent oscillations are 
confined in a plane defined by B and the z-axis. There- 
fore, the problem reduces to solving for the motion of 
the tilt angle (p of Q relative to the z-axis. Writing 
Q = Q (sin ip, 0, cos ip), we find 



ip = fiD, 

D = -ujQ sin(</> + 26»o ) 



(11) 



2. Inverted hierarchy with arbitrary fi 

Using the relation 9q = ir/2 — 9q, the potential |15 
can be written as 



V(<p) 



1 + cos((/3 - 29 ) 



h 
~2 



(16) 
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in the inverted hierarchy. Depending on the strength of 
the neutrino-neutrino interactions (i.e., the ratio /Lt/w), 
ip can take on a range of initial values from <p(0) ~ to 
28o — 7r in the small 9q limit. In other words, the evolution 
of the system begins with ip sitting near the maximum of 
the potential V(tp) if n/u> > 1, or, for smaller values of 
/i/cj, somewhere further down the slope. 

Since c^j(O) = 0, the form of the potential V(ip) guar- 
antees that ip always rolls to the minimum p m - ln ~ —it in 
the small 9q limit. Evaluating for P z and P z at p — tpmim 
one finds 



Z I <£n 



PA 



1, 



1, 



uj < 2/i, 

10 > 2/i. 



(17) 



Thus complete flavour conversion, i.e., P 2 = P z = — 1, 
is only possible in the strong neutrino-neutrino coupling 
limit ji/u) 1. A partial conversion can be achieved for 
comparable /i and uj. No conversion occurs for u) > 2/i, 
which remains always in the small oscillations regime. 



3. Inverted hierarchy with fi > u 

We now focus on fx/uj > 1. The initial condition 
tp(0) rs —(uj/hQ) 29q puts near the maximum of the 
potential V(y>) at < = 0. Since <p(Q) — 0, and <p(0) and 
26*o are both small, the motion of <p begins slowly, which 
explains the long plateau phases between dips in Fig. [Tj 
The dips themselves correspond to the crossing of the 
anharmonic potential once ip grows to the order unity. 
The duration of the dips is fixed by the crossing time of 
the anharmonic potential and thus is of order kT 1 . The 
duration of the plateau, on the other hand, is determined 
by the smallness of the mixing angle 9q. If 9q is exactly 
zero, i.e., no mixing, then there is no motion and p sits 
at the exact maximum of the potential forever. 

In order to estimate the time it takes for the polarisa- 
tion vectors to flip in the small 9q limit, we return to the 
equation of motion for this case, 



p = k 2 sin(p — 29q), 



(18) 



with 9q <C 1. As long as p is small this is equivalent to 
p = K 2 (p- 29 ). (19) 

Using the initial conditions <^(0) = — (oj/fiQ) 29q and 
ip(0) = 0, this equation is solved by 



tp(t) = 29 



1 - 1 + — cosh(Ki) 
HQ 



(20) 



Initially p(t) - tp(0) = -6~ (1 + uj/ fiQ)(nt) 2 , but at t 
of order n^ 1 turns to exponential growth when tp has 
become of order 9q. Therefore, the time it takes for p 
to grow to order unity, or equivalently, the half period of 
the bipolar motion is 



^bipola 



-K^ln^o (1 + u/tiQ)]. 



(21) 



Therefore, the duration of the plateau phases in Fig. [T] 
or the time between dips scales logarithmically with the 
small vacuum mixing angle. 



E. Which mass hierarchy? 

We demonstrated that, assuming small mixing, com- 
plete flavour conversion occurs for the inverted mass hier- 
archy, while small oscillations occur for the normal hier- 
archy. However, this applies only if the initial ensemble 
consists of v e and v e . If the initial ensemble had con- 
sisted instead of and then the situation would be 
reversed so that large flavour conversions would occur 
for the normal hierarchy. Put another way, the unstable 
case is when the initial ensemble consists of that flavour 
which is dominated by the heavier mass eigenstate. This 
symmetry between the neutrino flavours is an important 
difference to flavour conversion caused by the MSW effect 
that occurs for the normal mass hierarchy independently 
of the flavour of the initial state. 



F. Which flavour conversion? 

We have used the term "flavour conversion" loosely to 
describe the simultaneous conversion of equal numbers 
of v e and v e to equal numbers of and v^. Of course, 
the net flavour lepton number of the initial state vanishes 
and remains so at all times. Therefore, the "conversion" 
we are considering in the bipolar context does not violate 
flavour lepton number, but, rather, we should think of it 
as a coherent pair process of the form v e v e — > v^v^. 

If the initial system is asymmetric with a net electron 
lepton number, i.e., P and P are not initially identical, 
this difference is quantified by our vector D, whose pro- 
jection in the z-direction represents the net flavour lepton 
number. From Eq. ((7j) we observe that D-B is conserved. 
In other words, there is no net conversion of flavour be- 
yond what is caused by ordinary vacuum oscillations, and 
the net lepton numbers of vacuum mass eigenstates are 
strictly conserved. This statement is independent of the 
strength of /i, and applies irrespective of the asymmet- 
ric system being in the synchronised or bipolar regime. 
It also applies to the multi-mode system described in a 
later section. 



III. BACKGROUND MATTER 

An ordinary background medium has little impact on 
the bipolar flavour conversion [H, [23|, |24|. This sur- 
prising observation runs against the intuition that in a 
medium the mixing angle should be suppressed. We have 
already found that the time scale for bipolar flavour con- 
version depends only logarithmically on the vacuum mix- 
ing angle. One would perhaps expect this time scale to 
depend also logarithmically on the matter density. 
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To investigate this case we include matter effects 
caused by charged leptons in the equations of motion, 



d t P = [+loB + XL + n (P - P)] xP ; 
d t P = [-wB + AL + /i (P — P)] x P . 



(22) 



In the absence of other terms, P and P would pre- 
cess around L in the same direction with a frequency A. 
Therefore, it was noted in Refs. [13, HH that the equa- 
tions simplify if we study them in a frame co-rotating 
around the L-direction, i.e., around the ^-direction. 

In the co-rotating frame, the equations of motion for 
this system take the form of our original equations |T]), 
except that B is now time dependent, rotating around 
the z-direction with frequency — A, 



/sin(26» )cos(-At)~ 
B = sin(20 o )sin(-Ai) 
V - cos(20 o ) , 



(23) 



If this rotation is faster than all other frequencies one 
would naively expect that the transverse components of 
B average to zero, leaving us with (B) along the z-axis, 
i.e., an effectively vanishing mixing angle and no flavour 
conversion. However, as it turns out, there remains a 
net effect on the polarisation vectors, and bipolar flavour 
conversions occur after all! 

To understand this case quantitatively and qualita- 
tively we consider the full set of equations ([I]), keeping 
in mind the time dependence of the B vector in the co- 
rotating frame [cf. Eq. (|2"3"|) ]. We define a new vector Q 
as per Eq. ([6]). However, since B is now time-dependent, 
the new equations of motion have a slightly more complex 
structure than Eq. (0, 



Q 



(iDxQ B, 

A* 



D = wBxQ. 



with 



B = Asin(20 o ) 



sin(— Xt) 
- cos(— At) 




(24) 



(25) 



Thus Q is not strictly conserved; its length |Q| will, in 
the small 0o limit, exhibit order 9q fluctuations around 
some mean value given by Eq. ([5]). 



Note that D • B too is no longer exactly conserved. 
However, this non-conservation is likely to become rel- 
evant only in the vacuum oscillation dominated regime 
and for significant vacuum mixing angles: Since the non- 
constant part of B is proportional to sin(20o) and oscil- 
lates with frequency A, while, according to Eq. (|2"4")1 , D 
evolves with frequency w, the time varying part of D • B 
is proportional to sin(20o) and will tend to average out 
in the matter dominated regime we are interested in. Of 
course, if the matter term varies adiabatically, a large 
change of D • B in the form of the ordinary MSW effect 
is possible. 

We are interested in the case of small mixing in an 
inverted hierarchy, so that 



B 



20 o cos(-Ai) 
20 o sin(-At) 



(26) 



to the lowest order in 9q . Furthermore, since we are con- 
cerned only with instances at which the deviation of Q 
from the z-direction is small, we can parameterise the 
motion of Q by two small tilt angles tp x and tp y , so that 
to lowest order 



(27) 



and |Q| vanishes. As a consequence of this expansion, 
D z and D z are of higher order. In other words, D to 
lowest order does not develop a z-component; it remains 
a vector in the x-y-plane. One then finds a simple set of 
equations of motion, 




if x = K 



ip x - 20 o 1 



A 2 



<Py + 20q 1 



A* 2 Q 2 



A* 2 <9 2 



cos(Ai) 



sin(Ai) 



(28) 



with n 2 = ujuQ. Of course, had we considered the normal 
hierarchy in which B z w —1, we would have found two 
driven harmonic oscillator equations instead. 

Using the initial conditions (f x (0) = —(uj/hQ) 29q, 
^(0) = 0, p„(0) = 0, and p„(0) = (Aw/mQ) 20 o , 
Eq. (J2SK is solved by 



<p x (t) = -20c 



<p v (t) = -20 o 



K 



K 2 + X 2 
2 



1 + — ) cosh(Kt) - ( 1 



M 2 Q 



cos(Xt) 



K 2 + X 2 



- ( 1 + -^-] sinh(Kt) + ( 1 - -^-) sm(Xt) 



(29) 



Therefore, the tilt angles have a small oscillatory mo- tion driven by the rotating B vector. For t > k 1 these 
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oscillatory terms no longer matter much relative to the 
exponentially growing terms which scale asymptotically 
as e Kt . The most remarkable feature, however, is that 
the exponential term for (p y involves an additional factor 
X/k that is absent for tp x . 

If the matter effect is small, A <C k = (ui/iQ) 1 ^ 2 , the tilt 
is mostly in the x-direction in the co-rotating frame. The 
tilt in the y-direction is relatively suppressed by a factor 
X/k. In the limit A — > 0, the co-rotating frame coincides 
with the "laboratory" frame, and the tilt occurs exclu- 
sively in the x-direction as already seen in SecHU On the 
other hand, when the matter effect is strong, A>/s, the 
opposite applies. The tilt is mostly in the y-direction in 
the co-rotating frame, the motion in the ^-direction be- 
ing relatively suppressed by ft/A. We have observed this 
counter-intuitive behaviour also in numerical examples 
where indeed Q tilts in the y-direction when the matter 
effect is large. 

To track the flavour evolution, only the z-component 
of Q is of interest; its tilting direction is irrelevant. For 
t ^S> we ignore the oscillatory terms and use the 

asymptotic behaviour cos1i(k£) smh(nt) 



\e Kt so that 



vQt) = (<& + <£) 1 =0o 



(k 2 + A 2 ) 



K I U 



Therefore, the time scale for flavour conversion is 



^bipola 



-K _1 ln 



(k 2 + A 2 ) 



1/2 



(30) 



(31) 



The presence of matter has little impact on the overall 
behaviour of the bipolar system except for a logarithmic 
extension of Tbi po iar- This effect is illustrated in Fig. [21 
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FIG. 2: Evolution of P z and P z in several systems with 
background matter described by Eq. 1)220 , The parameters 
Go = 0.01 (i.e., inverted hierarchy), w = 1, and fj, = 10 
are common for all three systems. The blue/dotted line has 



A = 10 , the green/dashed line A = 
line A = 10 4 . 



10 , and the red/solid 



IV. VARYING NEUTRINO DENSITY 

In numerical simulations of the flavour evolution of su- 
pernova neutrinos in the single-angle approximation, one 
observes almost complete flavour conversion, apparently 
caused by the bipolar effect. We have seen that a bipo- 
lar system does lead to almost complete conversion, but 
also that the evolution is periodic. Therefore, being in 
the bipolar regime alone does not explain complete con- 
version. We have also seen that the impact of ordinary 
matter on the bipolar system is negligible. Therefore, 
it appears that the decline of the neutrino density, and 
therefore of /i, along the neutrino flux is the likely cause 
of almost complete flavor conversion. 

To study the impact of a time-varying (or, in a su- 
pernova, space-varying) neutrino density in a concrete 
example, we assume that all neutrinos have the same 
energy and oscillate in vacuum with the frequency u> — 
0.3 km -1 , corresponding to the "atmospheric" neutrino 
mass difference for typical supernova neutrino energies 
[see Eq. (|A9[) ]. We express all "frequencies" in units 
of km -1 and all length scales in km as appropriate for 
the supernova environment. For the v-v interaction en- 
ergy we use /x = 0.3 x 10 5 km -1 at the neutrino sphere 
(r = 10 km), and a dependence on the radius given by 
Eq. (|A10p . i.e., essentially an r~ 4 scaling. This scaling 
reflects both the ordinary flux dilution with r -2 , and the 
degree of collinearity in the neutrino motion which intro- 
duces approximately another factor r~ 2 . The quantity 
fi/u as a function of radius is shown in Fig. [3] 

Furthermore, we assume that (i) we have equal fluxes 
of neutrinos and antincutrinos, (ii) all of them are ini- 
tially in the same flavour, (iii) the mass hierarchy is in- 
verted, and (iv) the mixing angle is sin2#o = 0.001, rep- 
resenting a possibly small 0xa mixing angle. We then find 
numerically the survival probability of the initial flavour 
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FIG. 3: Neutrino-neutrino interaction strength /j, in units of 
uj for our toy model of supernova neutrino oscillations. 
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r (km) 

FIG. 4: Survival probabilities for v £ or v £ in our toy supernova 
model with symmetric initial conditions. 



as shown in Fig. |4j 

As expected, we observe in Fig. @] bipolar oscillations 
above the neutrino sphere at r = 10 km. Moreover, we 
observe that the oscillation amplitude declines as a func- 
tion of radius so that after a few tens of km we obtain 
complete flavour conversion. The question then is: can 
the decline of the upper envelope of the survival proba- 
bility be explained by way of the flavour pendulum? 

In the pendulum language, we start at small radii with 
the usual full oscillations. However, while the pendulum 
oscillates, we slowly reduce fx, i.e., we increase adiabat- 
ically the moment of inertia / = fx . Since the kinetic 
energy is D 2 jll and the angular momentum D is con- 
served, an increase in / corresponds to a decrease in the 
kinetic energy. This is akin to a dancer who can speed 
up or slow down a pirouette by changing the moment of 
inertia. The continuous decrease in /x and hence in the 
kinetic energy, however, means that during subsequent 
swings, the pendulum will be unable to reach its previous 
height. This explains the general feature of a declining 



upper envelope in the survival probability caused by a 
decreasing fj, (Fig.d]). 

It is important to note that the potential energy is in- 
dependent of /i; it depends only on u>. Therefore, our pen- 
dulum is not like a gravitational one where the potential 
energy would be affected by a changing mass. Instead, 
we can imagine the bob of the pendulum being charged 
and feeling the force of a homogeneous electric field. 

Let us now quantify the decline of the upper enve- 
lope in Fig. 3J The kinetic energy at a given time is 
T = fi(t)D(t) 2 /2. The conservation of angular momen- 
tum, except for the natural pendulum motion, implies 
that a sudden change A/i at some time t causes a change 
in kinetic energy of AT = Afj,D(t) 2 jl. For example, if we 
change /i by A/i when the pendulum swings past its low- 
est point at which the kinetic energy is maximal, then the 
relative change is AT max /T max = A/i//i. However, /i de- 
creases slowly compared to the oscillation period so that 
we may assume a linear decline. Therefore, the change 
of T max occurs over the entire oscillation period and thus 
must be weighted with a factor proportional to D{t) 2 
over one oscillation period. If the oscillation is approx- 
imately harmonic, we have D(t) 2 oc sin 2 (fti) with k the 
pendulum's natural frequency. The average of sin 2 (hit) is 
1/2. Therefore, if over one period fj, decreases by a factor 
(1 — e) where Kl, then AT max /T max = — e/2. In other 
words, T max is reduced by a factor (1 — e/2) = (1 — e) 1 / 2 
so that T max oc p}l 2 . 

The maximal kinetic energy equals the maximal po- 
tential energy minus its minimum, achieved within one 
oscillation. The potential energy normalised such that its 
minimum is zero, oj (2 + B • Q) for small mixing angles 
and for large fJ,/u>, gives us the projection of the summed 
polarisation vector S on the flavour axis. Therefore, the 
upper envelope of the survival probability in Fig.|3]should 
scale with /i 1 / 2 . In Fig. [5]we show the same survival prob- 
abilities with (/i/lu) 1 / 2 as a radial coordinate. The decline 
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FIG. 5: Survival probability for v e or v e as in Fig. [4] here 
plotted as a function of (/i/aj) 1,/2 . 
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is indeed nearly linear, especially for small-amplitude os- 
cillations (toward the right side of the plot) where the 
pendulum oscillations are more nearly harmonic. 

Our pendulum analogy elegantly explains the most 
puzzling feature of the bipolar supernova neutrino oscil- 
lations, i.e., that they actually lead to flavour conversion 
rather than permanent bipolar oscillations. We are sim- 
ply seeing the relaxation of the pendulum to its down- 
ward rest position as kinetic energy is extracted by the 
reduction of the neutrino-neutrino interaction potential 
and thus the increase of the pendulum's inertia. 



V. NEUTRINO ASYMMETRY 



A. Realistic supernova example 
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The behaviour of neutrinos streaming off a supernova 
core looks, however, rather different from this simple pic- 
ture. Besides the dependence of /x on the radius, another 
crucial feature is the initial neutrino-antineutrino asym- 
metry; the number flux of v e is the largest, that of v e 
smaller, and those of the other species yet smaller but 
equal to one another. We represent this situation by the 
initial conditions P z (0) = 1 and P 2 (0) = 0.8. The equa- 
tions of motion for this system are simply those given in 
Sec. Ill Al Solving them numerically yields the relative v e 
and v e fluxes shown in Fig. [5J Plotting relative fluxes in- 
stead of survival probabilities makes conservation of the 
net flavour-lepton flux evident. 

Observe that the initial flavour-lepton asymmetry is 
conserved so that there remains a net flux of v e origi- 
nally set at r = 10 km. Otherwise there is complete 
flavour conversion over a length scale given by the de- 
crease of [i as a function of radius. Similar results are 
found in detailed numerical studies within the "one-angle 
approximation" where the neutrino interaction strength 
is taken equal for all modes (Fig. 8c of Ref. [23| and pri- 
vate communication by S. Pastor and R. Tomas based on 
the numerical scheme of Ref. 29] ) . Changing the vacuum 
mixing angle and adding normal matter causes only the 
minor logarithmic changes predicted earlier. 

The behaviour of the neutrino ensemble between the 
neutrino sphere at 10 km and r k, 45 km is ex- 
plained by synchronised oscillations due to the neutrino- 
antineutrino asymmetry and the large value of /i in this 
region. Beyond this we enter the bipolar regime out to 
about 200 km when vacuum oscillations take over. In 
other words, the flavour evolution of neutrinos streaming 
off a supernova core are determined by a transition from 
synchronised oscillations at small radii (large n), bipolar 
oscillations at intermediate radii (intermediate fj,), and 
ordinary vacuum oscillations at large radii (small fi) as 
first stressed in Ref. (22[. If ordinary matter is included, 
it affects the synchronised region at small radii in the 
usual way by making the effective mixing angle smaller, 
and likewise at large radii where no collective effects oc- 
cur. It is only in the intermediate, bipolar oscillation 



FIG. 6: Relative fluxes of z/ e (blue/dotted) and i> e (red/solid) 
in our toy supernova model with 20% fewer antineutrinos than 
neutrinos and sin2#o = 0.001. 

regime where ordinary matter has no significant impact 
on the system. Confusingly, the bipolar behaviour does 
not correspond to the limit of large or small neutrino 
densities; it corresponds to intermediate densities [22|. 

We have already explained the decline of the upper en- 
velope of these curves in the bipolar regime which should 
scale as /x 1 / 2 . To confirm this behaviour once more we 
show in Fig. [7] (middle panel) the relative fluxes of Fig. [5] 
with (/i/w) 1 / 2 as a radial coordinate. The decline of 
both the upper and lower envelopes are stunningly lin- 
ear. This reflects the small-amplitude nature of the oscil- 
lations which are now nearly harmonic so that our pre- 
vious argument works better here than in the previous 
section. 



B. Transition between different oscillation modes 

The lepton asymmetry of the neutrino flux is crucial for 
understanding a realistic supernova case. To develop a 
first understanding of this situation we consider an asym- 
metric ensemble with the initial condition P(0) = aP(0) 
where < a < 1. If the v-v interaction is sufficiently 
strong, the two polarisation vectors will hang together 
by their "internal magnetic field" and all vectors P, P, 
S, and D precess around B with the same synchronised 
frequency. Using D = wB x S and D = S(l — a)/(l + a), 
we thus find the usual result 

1 + a 

^synch = ~ — W . (32) 

Observe that the synchronised motion is faster for more 
symmetric systems! 

To achieve synchronised oscillations, the internal fre- 
quencies of the system, i.e., their precession frequency 
around a common spin, must exceed the slow precession 



10 



w 

CD 
X 
3 



CD 

> 
• t—i 
+-> 

'CD 



1 

0.8 
0.6 
0.4 
0.2 

1 

0.8 
0.6 
0.4 
0.2 

1 

0.8 
0.6 
0.4 
0.2 






1 




10 5 



FIG. 7: Relative fluxes of u e (blue/dotted) and v e (red/solid) 
as in Fig. |SJ here using (p/cu) 1 ^ 2 as a radial coordinate. Top: 
Vacuum mixing angle sin26io = 0.1. Middle: sin2#o = 10~ 3 
as in Fig. [SJ Bottom: sin 26 = 10" 5 . 



around the external magnetic field w sync h. In our case of 
two polarisation vectors, the internal motion is described 
by S = /jD x S. When synchronisation prevails, D has 
a fixed length (1 — a) so that the internal frequency is 
(1 — a) fi. Successful synchronisation means that this 
frequency must exceed w syn ch, or, equivalently, 



1 



1 — a 



<(l-a)(i. 



(33) 



This is exactly the argument and result presented in 
Ref. (22J. When this condition is violated, the motion 
becomes bipolar. 

Actually, this condition derives from an even simpler 
argument if we consider the total energy of the system, 



H = wB • (P + P) 



P) 2 



(34) 



Synchronised oscillations require that the energy of the 
system be dominated by the spin-spin interaction, the 
second term. Assuming a small vacuum mixing angle 
so that B is nearly aligned with the initial polarisation 
vectors, and observing that |P + P| = l + a and |P — P| = 
1 — a, the requirement that the spin-spin term dominates 
is ui (1 + a) < (1 - a) 2 fi/2, identical with Eq. ([33]) up to 
a multiplicative factor. 

In other words, synchronised oscillations occur when 
the neutrino-neutrino part of the Hamiltonian always 
dominates over the vacuum-oscillation part, even for the 
most disadvantageous orientation of the polarisation vec- 
tors. At the other extreme, no collective effects obtain 
when the neutrino-neutrino part never dominates, even 
for the most advantageous orientation of these vectors, 
i.e., when \x < u>. The bipolar regime corresponds to 
the intermediate range where the relative magnitude of 
the different energy contributions depends on the orien- 
tation of the polarisation vectors. In summary, bipolar 
oscillations are expected when 



1 — aV 



(35) 



and thus occur for an intermediate strength of the 
neutrino- neutrino interaction, i.e., an intermediate range 
of neutrino densities as first discussed in Ref. [22}. The 
exact numerical factor on the r.h.s. of this equation is 
taken from Eq. (|46|) below. 



C. Asymmetric system as a pendulum with spin 

This simple reasoning gives us the correct scale for the 
transition between synchronised and bipolar oscillations, 
but does not explain the nature of the transition. The 
two oscillation modes must be extreme cases of a con- 
tinuum, yet the nature of the intermediate cases is not 
obvious. It turns out that this continuum has a straight- 
forward physical interpretation that is quite illuminating 
for an understanding of the entire system. 

To this end we note that the asymmetric system is de- 
scribed by the same equations of motion ([7]) for Q and D 
as the symmetric case, and Q = |Q| is likewise conserved. 
The new feature here is the initial condition D(0) 7^ 
and D • Q 7^ 0, so that Eq. ([7]) now implies 



D = - 



1 Q x Q D 

/' 



Q 



Q 



Q 



Q 



Since 



cr = D Q/Q 



(36) 



(37) 



is a constant of the motion, the expression (|36|) can be 
equivalently written as 



D 



q x q 



crq, 



(38) 
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with q = Q/Q. The first term in the expression corre- 
sponds to the orbital angular momentum as before, while 
the second term plays the role of an inner angular mo- 
mentum (i.e., spin) of the pendulum's bob. This spin 
is always along the direction q, implying that we should 
think of the system as a spinning top mounted in a way 
that its axis of rotation can swing like a pendulum. If 
the top has no spin (i.e., a — 0), it acts as an ordinary 
spherical pendulum. 

Starting with the equation of motion D = wB x Q and 
observing that q x q = 0, we find 

— + ^q = wQBxq. (39) 
M 

For the case a = 0, we recover the equation of motion 
of a pendulum swinging in the plane defined by B and 
the initial vector q(0) which we always choose to be in 
the z-direction. Conversely, if [i is so large that we can 
neglect the first term, we are back to a spin-precession 
equation with a precession frequency w P rcccss = ojQ/ct. 
For /i> uj we have 

Q « 1 + a, 

a w I- a, (40) 
so Wprocess is indeed equal to w S ynch 

of Eq. (|32|). 

The extreme cases thus have an intuitive interpreta- 
tion. Synchronised oscillations correspond to the flavour 
top spinning so fast that its response to the force field is 
precession, just like a spun-up top on a flat table surface. 
On the other hand, if the top spins slowly (correspond- 
ing to the case of a small neutrino asymmetry), then it 
swings like an ordinary pendulum; the inner angular mo- 
mentum in this case has little impact, and we are in the 
bipolar mode. 

Let us consider the asymmetric system in more detail. 
Several conserved quantities are apparent. One is the 
energy of the system, 

E = £7pot + -Ekin 

= ^Q(B-q+l) + |D 2 

= wQ(B-q+l) + ^q 2 + !<7 2 , (41) 

where we have added a constant to the potential energy 
such that it vanishes when the pendulum is oriented op- 
posite to the force field. The other conserved quantity is 
the projection of D in the B direction, corresponding to 
the conservation of that component of the angular mo- 
mentum parallel to the force field and which is thus not 
subject to a torque. 

The conservation of angular momentum implies that 
the initial lepton asymmetry cannot be changed beyond 
the amount caused by ordinary vacuum oscillations. In 
all practical cases we begin with D oriented along the 
z-axis, while B is tilted by 20$. If 6q is small, the initial 
neutrino asymmetry is almost perfectly conserved. Thus 



the self-interacting system cannot stimulate an exotically 
large flavour conversion effect. 

If the spin is not quite fast enough for perfect pre- 
cession, the overall motion is a wobble. The spinning 
top starts in a nearly upright position, its axis pointing 
nearly to the north pole. It will tilt a bit until it reaches 
a certain latitude, when its motion reverses back to the 
north pole. The latitude of reversal will lie further south 
if either or both of a and /i is smaller. In other words, 
a smaller v-v term makes the system "more bipolar" . A 
more symmetric system (smaller a) is also more bipolar. 

The southernmost position the pendulum can reach is 
defined by energy and angular momentum conservation. 
As an example, we take the mixing angle to be very small 
so that B is very close to the z-direction. In this approx- 
imation, the initial total energy is E — 2ojQ + [i<j 2 /2 and 
energy conservation implies 

q 2 

wQ(l-cosw) = — . (42) 
2/i 

On the other hand, angular momentum conservation 
along the B-direction gives 

a = a cos if + /i -1 ^ sin ip , (43) 

where q± is the velocity perpendicular to B, i.e., the 
pendulum's velocity along a circle of latitude. The largest 
excursion </? max is reached when the pendulum reverses 
its motion at its southernmost position where q 2 = q\. 
Combined with Eqs. g2]) and (|33"]l. we find 

ua 

cos tp max = — — - 1 (44) 
2loQ 

for the largest excursion angle. 

As expected, if either or both of fi and a becomes 
smaller, Eq. (f44|) tells us that the pendulum reaches more 
southern latitudes. On the other hand, the equation has 
no solution for 

[ia 1 > AluQ . (45) 

This condition corresponds to the fully synchronised case, 
which prohibits any deviation from perfect z-alignment 
because of our artificial assumption of B and the initial 
P and P being exactly aligned. If the pendulum had not 
initially been perfectly aligned with B, solutions would 
exist for all values of the parameters. Still, for small 
mixing angles, Eq. (|45|) provides an excellent estimate of 
the condition for synchronised behaviour. Using Eq. (I40[) . 
this condition is equivalent, in the (j,/uj 3> 1 limit, to 

(!-a) 2 „w 

\ ~ >4~ (46) 
1 + a ^ 

where, we recall, a parameterises the lepton asymmetry 
by virtue of |P| = a|P|. 

Taking our previous asymmetric example with P z (0) = 
1 and P z (0) = 0.8 (i.e., a = 0.8) and assuming fj, ^> u, 
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synchronised behaviour is expected for fj,/u> > 180 or 
(/i/cj) 1 / 2 > 13.4. This estimate corresponds well with the 
onset of synchronisation in the top panel of Fig. [JJ where 
the mixing angle is large. Of course, the true point of 
onset also depends logarithmically on the mixing angle — 
see the other panels of Fig. [7J 

It is now evident that the onset of bipolar oscillations 
does not imply full conversions. As we move into the 
bipolar regime, the spinning top begins to wobble, reach- 
ing only some southern latitude, but not the south pole. 
How far south it will get, i.e., the spread between the 
upper and lower envelopes in Fig. [7J depends on the de- 
tails of how the system enters the bipolar regime. If the 
mixing angle is large, bipolar oscillations begin almost 
immediately so that the amplitude of the oscillations will 
be small. If the mixing angle is small, the delayed on- 
set of the first bipolar swing allows /i to decrease fur- 
ther, thereby resulting in a more southern turning point. 
Therefore, smaller mixing angles imply a later onset of 
oscillations and a larger spread between the envelopes. 
This is borne out by the examples shown in Fig.[jj where 
sin 2^o = 0.1, 10~ 3 and 1CT 5 from top to bottom. 



D. Equipartition of Energies 

We have noted that the energy of the spinning top de- 
creases in proportion to y}l 2 once it has entered the bipo- 
lar regime, assuming the decline of /i is sufficiently slow. 
It is illuminating to note that from that time onward, the 
total energy Eq. (|4Tj) is equipartitioned between E pot , the 
potential energy in the external force field, and -E^in, the 
internal and orbital kinetic energy of the spinning top 
(equivalent to the neutrino-neutrino interaction energy). 

To illustrate this point we show for our toy supernova 
the evolution of E pot , E kin and E tot /2 = (E pot + E kia )/2 
in Fig. [51 using once more (/x/w) 1 / 2 as a radial coordi- 
nate. Observe that indeed E pot = -Skin at /j,/uj = 180 or 
(/i/cl)) 1 / 2 = 13.4 for our example a — 0.8. It is intrigu- 
ing that the transition between the synchronised and the 
bipolar regime is practically independent of the mixing 
angle. The same is true also for the total energy E to t, 
which decreases very nearly as /it 1 / 2 in the bipolar regime 
as explained earlier. In the synchronised regime close to 
the supernova, the potential energy does not depend on /i 
whereas the kinetic energy decreases with [i. In the bipo- 
lar region, on the other hand, the kinetic and potential 
energies are nearly equipartitioned after averaging over 
the pendulum's nutation period. 

To understand this equipartition effect analytically, we 
consider a case where the nutation amplitude is very 
small as in the top panel of Fig. [8J (large mixing angle) , 
so that it suffices to study only the precession, i.e., we 
assume the pendulum's orbital motion is such that its 
velocity is along a circle of latitude. As a function of 
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FIG. 8: Different energy components E po t, -Ekin and Etot/2 
for the schematic supernova model of Fig. [U using (/i/aj) 1 / 2 
as a radial coordinate. The panels from top to bottom have 
the vacuum mixing angle sin2^o = 0.1, 10~ 3 and 1CP 5 , re- 
spectively. The vertical line at fi/u = V 180 ~ 13.4 marks the 
transition to the bipolar regime according to Eq. (|46[1 for this 
example where a = 0.8. 

excursion angle tp, the total energy is 

E tQt = ujQ{\ + cos tp) + |- + ^a 2 . (47) 

Using angular-momentum conservation Eq. (|43[) and the 
relation q 2 = q\ to eliminate the orbital velocity, we find 

E tot =uQ(l + cos ip) + ^- • (48) 

2 1 + cos ip 

We further recall that the system enters the bipolar 
regime when /ia 2 = 4u)Q, and that at this point the ex- 
cursion angle is still small so that cos ip = 1. Therefore, 



13 



Spot = £"kin = 2uiQ at the onset of the bipolar regime. 
Subsequently, E tot is expected to scale as fj, 1 ^ 2 so that 
Etot = {^ojQ^ia 2 ) 1 / 2 . We can now solve for the expres- 
sion (1 + cost/?) as a function of fi and find explicitly 



E, 



pot 



Skin = (luQh<7 2 ) 



2\l/2 



An important detail is that equipartition cannot hold 
all the way to very small /i. Angular momentum con- 
servation, i.e., the approximate conservation of the net 
is e flux in the limit of a small vacuum mixing angle, im- 
plies that the potential energy is bounded from below. In 
terms of polarisation vectors, this means that the strict 
conservation of B • (P — P) leads to an approximately 
constant P z — P z = 1 — a = a in the case of small vac- 
uum mixing. Therefore, the smallest allowed value of the 
potential energy is 



Spot > 2(1 — a)uj = 2auj . 



(49) 



In the example of Fig. [5] we have used an asymmetry a — 
0.8. This gives an analytic estimate of E pot > 0.4 w using 
Eq. (|49|) . which is in excellent agreement with numerical 
results. On the other hand, the kinetic energy, being 
proportional to [A, must eventually vanish. Therefore, 
-Epot and Ek nl approach different limits as (i — > 0, as 
borne out by Fig. [8] 

Therefore, the evolution of the system now appears in 
a different light. When /i decreases slowly, the system 
never properly enters the bipolar regime: it stays at the 
edge of it. The polarisation vectors spiral and slowly 
tilt in such a way that their energy in the external B- 
field and the internal spin-spin energies stay equal to the 
degree allowed by net lepton number conservation. 



VI. MANY MODES 
A. Multiple frequencies 

We now turn to a more realistic case of an ensemble of 
v e and v e with different energies, i.e., different vacuum 
oscillation frequencies w,. The equations of motion are 

dtPi = [+WiB + (i (P - P)] x P t , 

d t P l = [-UiB + n (P - P)] x P, , (50) 

where for N modes we use 

JV JV 

p=J2 p i and p = J2 Pi ( 51 ) 

i=l i=l 

We keep the normalisation |P| = |P| = 1 for the entire 
ensemble so that the individual modes are normalised to 
1^1 = 1^1=^-1. 

In full analogy to the previous treatment we introduce 
the vectors Si = P, + Pi and D; = Pi — Pi as well as 
D = ^ Dj = P — P so that 



Si = LUiB xD, + £iD x Si , 
Di = w-iB x Si + x Di . 



Each pair of modes Pi and Pi evolves symmetrically so 
that each Di is always oriented along the y-axis and the 
terms D x Di vanish. 

We now assume strong coupling with [ijuji 3> 1 for all 
modes so that we can also drop the WiB x Di term. This 
leaves us with the approximate equations of motion 



Si = /iD x Si 
D, = u>,B x S,: 



(53) 



Equation (|53|) implies that all Si evolve in the same way 
because they precess in the same field D. Flavour conver- 
sion is now described by a single vector S = Si (and 
thus Si = S /N) , and governed by 



S = fiD x S , 



(54) 



Thus, the evolution of the flavour content proceeds in 
the same way as before [cf. Eq. ©], but with the role of 
lj replaced with the average oscillation frequency of all 
modes (ui) = N^ 1 uii. 

On the level of the individual modes, the "tilting mo- 
tion" around the y-axis is the same for all neutrino and 
antineutrino modes, so in this sense their motion is syn- 
chronised. On the other hand, the transverse motion 
characterised by Di is different for every mode because 
before summing, the equations of motion are 



S 
D, 



/iDxS, 
B x S 



N 



(55) 



(52) 



where we have used Si = S/N. 

It is, therefore, incorrect to say that all modes form two 
block spins P and P which evolve separately in the bipo- 
lar sense, with each individual mode staying aligned with 
its respective block spin. Separate alignment for neutri- 
nos and antineutrinos was explicitly claimed, for exam- 
ple, above Eq. (4.1) in Ref. [12|, in reference to the au- 
thors' own numerical studies of Refs. (hIEI- However, 
the observation of separate alignment does not appear to 
be documented or demonstrated in these papers. What- 
ever the origin of these authors' observation of bipolar 
alignment, it is in conflict with our analytic treatment. 
We have numerically verified that bipolar alignment does 
not hold and that our individual Pi and Pi vectors do 
indeed evolve differently, as shown in Fig. [51 

The z-components of all modes evolve identically while 
the transverse motion is different for modes with different 
LUi. The transverse motion for the neutrino and antineu- 
trino polarisation vectors are opposite so that neutrinos 
and antineutrinos form two distinct cohorts. In this sense 
the evolution actually is bipolar. Therefore, we stick to 
this established terminology, keeping in mind a broad in- 
terpretation of the word "bipolar." 

The most important conclusion is that the strongly in- 
teracting multi-mode system is exactly equivalent to one 
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FIG. 9: Evolution of the individual modes P y (red/solid) and 
P y (blue/dotted) in a 4- mode system described by Eq. (I5U|I . 
with 6(i — 0.01 (i.e., inverted hierarchy), = 1,2,3,4 for 
i — 1, 2, 3, and 4, respectively, and strong neutrino-neutrino 
interaction fi — 100. 



mode if one interprets the vacuum oscillation frequency 
lj as the average of all modes. Therefore, the entire sys- 
tem is still characterised by a single collective variable. 
This conclusion also holds for the asymmetric system of 
unequal densities of v and v. The variation of different 
vacuum oscillation frequencies is not a source of kincmat- 
ical decoherences and the system behaves, in this sense, 
similarly to synchronised oscillations. 



B. Different interaction strengths 

Instead of different frequencies one may also consider 
different coupling constants between different modes. In 
this case, the general equations of motion are 



d t V,, = 



d t V,, = 



+cJiB+J2vij (Pi -Pi) 

3=1 
N 

~,Ti • ^ ,/,,(!», I\) 



3=1 



X P, 



xPi. (56) 



Variations in the interaction coefficients are motivated 
by their dependence on the relative angle of the neutrino 
trajectories. Neutrinos streaming off a supernova core 
are far from isotropic so that unequal /ijj coefficients are 
unavoidable. The recent multi-angle simulations were 
aiming precisely at this issue [23l . |24| . 

Even in an isotropic medium the coupling constants 
between the different modes are never equal, but involve 
a factor (1 — cos# pq ) from the current-current structure 
of the weak-interaction Hamiltonian. However, isotropy 
implies that all modes with different momenta p but iden- 



tical p = |p| evolve in the same way. Therefore, the an- 
gular part of the integral in Eq. (|A2[) can be trivially 
performed in the sense that the cos # pq term averages to 
zero. The isotropic system is thus equivalent to the case 
of multiple frequencies, but a common coupling constant 
described in Sec. IVI Al In other words it is equivalent to 
the "single-angle" case. 

In terms of our variables S; (the flavour-dependent par- 
ticle plus antiparticle number) and Di (the net lepton 
number), the equations of motion are 



Si = ~,B xD, + ^ /i,J), x Si , 

3 = 1 
JV 

I), = ~ ,B X Si + ^2 Mi3 D 3 x D * ■ 



(57) 



3 = 1 



For the global S and D vectors this implies 



N 



N 



S = ^cjjB x Dj + ^2 A% D 3 x Si , 



'■3 = 1 



N 



D = ^WjB x Si 



(58) 



where we have used the symmetry fj^j = fj,ji. Equation 
(fS"8")) cannot be brought into the form of a closed set of 
equations. However, even in this general case the quan- 
tity B • D is a constant of the motion. This means that 
bipolar oscillations never lead to lepton number conver- 
sions beyond the amount caused by vacuum mixing. 

It is possible to formulate a sufficient condition for a 
multiple-coupling constant system to behave as a simple 
flavour pendulum. Starting with Eq. (|5T[) we note that 
the system acts as a single flavour pendulum if all Si 
and Di tilt with the same speed. Self-consistency then 
requires that all frequencies are equal, u>i = uj, and that 



N 

E 



fJrij — M : 



(59) 



with fi a number that is independent of i. In addition, 
symmetry between different modes requires that /Lty = 
fiji. These conditions are met, for example, if 



Hij oc g 



N 



(60) 



where g(x) is an even function that is periodic in the 
sense g(x + 1) = g[x). An example is g{x) = cos(27ra;). 

An important case where these conditions are violated 
is neutrinos radiated from a supernova core. Assuming 
overall spherical symmetry, the only parameter that dif- 
ferentiates between trajectories is cos 9 with 9 the angle 
relative to the radial direction. Considering instead the 
schematic case of neutrinos emitted "isotropically" by a 
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plane surface, the coupling constants are [23|], following 

Eq. Rxm . 



^3 



1 - 



N 2 



(61) 



Here, Cj = cos #i is the cosine of the neutrino mode rel- 
ative to the radiating surface's normal direction. We as- 
sume a uniform distribution in < c$ < 1 represented 
by discrete modes as in the second line of Eq. (JBXJ) • Note 
that ^i=i — |) = N 2 /2 so that the average (fcj) = [i 
is exact, even for a small number of modes. 

This example does not satisfy the condition Eq. (|5T 
because 



N 



4A 



1 - 



2N 



(62) 



(Note that the overall factor of N is compensated by us- 
ing individual polarisation vectors that are normalised to 
the length TV -1 .) Therefore, one cannot expect neutrinos 
streaming off a supernova core to oscillate in a collective 
manner. Rather, one should expect kinematical decoher- 
ence within a few bipolar periods, i.e., on a time scale 
of order kT 1 = (uj/j) 1 / 2 . In other words, the length of 
the total P or P vector is no longer conserved and the 
ensemble partly or fully decoheres. 

In simple numerical examples of a symmetric system 
this expectation is indeed borne out, i.e., a non-isotropic 
ensemble consisting of equal numbers of neutrinos and 
antincutrinos turns into an equal mixture of both flavours 
within a few bipolar oscillation periods for both the nor- 
mal and inverted hierarchy. 

On the other hand, the large-scale numerical studies of 
Ref. [HI show that neutrinos streaming off a supernova 
core are sometimes quite well represented by the single- 
angle case, i.e., collective behaviour rather than quick 
kinematical decoherence appears to be more generic. 
Again, the main differences to a simple symmetric system 
are twofold: there is a neutrino-antineutrino asymmetry 
and the effective density declines with radius. We have 
already observed in Sec. IV Dl that for the inverted hier- 
archy the asymmetry, together with the slow decline of 
the neutrino density, has the effect of slowly turning the 
polarisation vectors without the system ever entering the 
bipolar regime, i.e., the system teeters along the edge of 
the bipolar condition. Since the bipolar regime is never 
properly entered, it is less surprising that kinematical 
decoherence is not a prominent feature of the evolution. 

A dedicated research project is required to develop a 
deeper understanding of the conditions that determine if 
the system evolves as as single collective system in the 
form of the flavour pendulum, or if it kinematically de- 
coheres in that the individual polarisation vectors of the 
different modes are "randomised" in flavour space. 

Our main conclusion is that different oscillation fre- 
quencies are not a source of kinematical decoherence, 



while the multi-angle nature of a non-isotropic system is 
such a source, especially for a symmetric system. The de- 
tailed interplay between the collective mode represented 
by our pendulum and the multi-modal nature of the non- 
isotropic system remains to be understood. 



VII. CONNECTION TO QUANTUM PHYSICS 
A. Quantisation of the flavour pendulum 

We have seen that the flavour conversion of neutrinos 
streaming off a supernova core can be understood almost 
completely when we cast the equations of motion in the 
form of a pendulum in flavour space that may include 
inner angular momentum if we need to account for a lep- 
ton asymmetry. One may then ask if flavour conversion 
was possible in the absence of flavour mixing since, even 
if the pendulum were placed exactly on its tip, Heisen- 
berg's uncertainty relation would prevent it from staying 
there forever, just as an idealised pencil cannot stand on 
its tip indefinitely. 

To estimate the relevant time scale we recall that the 
equations of motion for the pendulum's excursion angle 
can be derived from the classical Hamiltonian Eq. (fT4]) 
with ip and D as the canonical variables. In order to 
quantise this system, however, we need to be more careful 
about absolute scales, and the equations of motion for the 
quantum variables must follow from the Hamiltonian of 
the full quantum system. In the case of N neutrinos and 
N antincutrinos, the full Hamiltonian is simply N times 
the one of Eq. (fTI)) , 

N Nk 2 

H tot = —fxD 2 + [l-oosr> + 20 o )l. (63) 

2 [i 

Identifying ip as the canonical coordinate, the familiar 
equations of motion (fTTj) follow classically using Hamil- 
ton's equations, provided we interpret ND as the con- 
jugate momentum. The pendulum's potential energy in 
the macroscopic sense scales with N, and likewise its mo- 
ment of inertia / = N/fi. On the quantum level, the 
corresponding commutation relation is 



N 



(64) 



Using this commutation relation, the same equations of 
motion (|11[) follow quantum mechanically from the full 
Hamiltonian (|63|) by way of Heisenberg's equations of 
motion ihip = [ip, H to t] and ihD = [D,Htot\- 

The fact that the same equations of motion for the 
tilt angle tp follow both classically and quantum mechan- 
ically irrespective of the size of N — provided we identify 
the appropriate canonical momentum — indicates that the 
flavour evolution does not depend on the size of the sys- 
tem. Thus, as long as our calculation is classical, we 
can work with polarisation vectors and associated angu- 
lar momenta that are normalised to unity. On the quan- 



16 



turn level, however, the absolute length of the polarisa- 
tion vectors will affect the quantisation of the system, 
and it is necessary that we use the correct Hamiltonian 
with the appropriate factors of N. 

To estimate the time scale for the pendulum to stay 
upright we consider its downward vertical position, i.e., 
we consider it to be a harmonic oscillator in its quantum- 
mechanical ground state. The uncertainties of the canon- 
ical variables in this state are 

{ip2) = \TK> and ^ ND ) 2 ) = \ 1 ^ ( 65 ) 

where n is the oscillation frequency of the pendulum and 
I = N / n its moment of inertia. In the strong-interaction 
limit fi 3> to, the expression for (if 2 ) becomes 

Let us now put in some realistic numbers. The typi- 
cal density of neutrinos in the supernova region of inter- 
est is estimated to be 10 32 cm~ 3 in Appendix [A] before 
Eq. (|A10[) . The volume of the critical region may be 
of order 100 km 3 so that some 10 53 particles may be 
in the system at any one time. Moreover, the bipolar 
regime begins at (/i/w) 1 / 2 of order 10. Thus a typical 
value for (t^ 2 ) 1 / 2 would be of order 10~ 26 . If this num- 
ber is taken to be a typical excursion angle caused by 
quantum fluctuations, then the time scale to tilt will be 
of order K _1 ln(10 26 ) » 60 n . Since k _1 is a frac- 
tion of a km, the quantum effect would happen over a 
length scale of tens of km. Of course, the vacuum mix- 
ing angle and thus the initial excursion of the pendulum 
is much larger than this quantum estimate. Moreover, 
the system would be subject to other forces. Still, this 
estimate demonstrates that in an unstable system expo- 
nential growth can quickly enhance quantum effects to a 
macroscopic scale. 

B. Full quantum system 

The discussion above shows that our classical treat- 
ment of the flavour evolution of a large neutrino ensem- 
ble was justified. However, it is still instructive to briefly 
explain the structure of the full quantum Hamiltonian. 
It is well known that the neutrino interaction in flavour 
space has an SU (2) structure and as such is equivalent 
to a spin system [3l], HH, [33[ . The equations of motion 
with a matter term ([!?!?]) can be shown to follow from the 
quantum Hamiltonian 

H = H + iJmattcr + H vv (67) 

= luB • (S - S) + AL • (S + S) + ii(S + S) 2 , 

where we now use carets to denote quantum operators 
explicitly. The operators S and S each represent an an- 
gular momentum N/2, i.e., N spin 4 particles that are 



linked to form one big "block spin" each. The equation of 
motion for S follows from ihdtS — [S,H] and the angu- 
lar momentum commutation relation [Si,Sj] = iheijkSk 

and similarly for S. Note that % drops out of the spin- 
precession equation dtS = wB x S. Therefore, spin pre- 
cession is fundamentally a classical phenomenon and one 
does not need to distinguish carefully between equations 
of motion for quantum operators and for expectation val- 
ues. However, for the nonlinear neutrino-neutrino term, 
it is not intuitively obvious that one can ignore correla- 
tion effects when taking the expectation values [HI, [H[ . 

The connection to the polarisation vectors is that P = 
(S) (2/iV) is the normalised expectation value of the spin 

that represents the particles, whereas P = — (S) (2/-/V) 
includes a minus sign. In other words, the quantities P 
and P play the role of "magnetic moments" in flavour 

space, whereas the quantities S and S play the role of 
angular momenta. We call them "flavour spins," but 
the terminology "neutrino flavour isospins (NFIS)" has 
also been used [IH, H3, [24| . The negative sign between 
the polarisation vector and flavour spin for antineutri- 
nos is consistent with antiparticles of equal spin carry- 
ing negative magnetic moments relative to the particles, 
such as the case of electrons and positrons. This nega- 
tive sign also explains that under the mass Hamiltonian 
in vacuum Hq, neutrinos and antineutrinos precess "in 
opposite directions" in flavour space. We find that the 
language of polarisation vectors is useful in the classi- 
cal limit, whereas the language of flavour spins is useful 
when dealing with the quantum aspects of the system. 

We note that only the vacuum Hamiltonian Ho distin- 
guishes between neutrinos and antineutrinos, while the 
matter and v-v parts of the full Hamiltonian can be ex- 
pressed in terms of a single big angular momentum oper- 
ator J = S + S. Therefore, these parts of the Hamiltonian 
are equivalent for our system consisting of neutrinos and 
antineutrinos, and a neutrino-only system consisting of 
two flavours. In the neutrino-only case, "flavour spin up" 
means v e , "flavour spin down" v^. In our case, "flavour 
spin up" means either v e or P M , and "flavour spin down" 
v e or , and the states v e and are fully equivalent in 
the absence of Hq (and similarly for v e and v^). 

In terms of the big angular momentum operator J, we 
see that the neutrino-neutrino Hamiltonian is of the form 
J 2 , while the interaction with ordinary matter is propor- 
tional to J z . These two operators commute so that they 
have a common set of energy eigenstates. This observa- 
tion is the quantum analogue to our classical result that 
the presence of ordinary matter leaves bipolar oscillations 
nearly unaffected. 

Some time ago it was speculated that a system of many 
spins interacting by a nonlinear Hamiltonian of the form 
J 2 could exhibit quantum entanglement effects in the 
sense that its evolution is coherently accelerated [35j j . 
Applied to our case, this conjecture means the follow- 
ing. Consider a "dense gas" consisting of exactly one 
v e and one v e with the same density as our macroscopic 
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system, i.e., with the same spin-spin interaction energy 
/i. In this case the four possible states of the system 
are grouped into a triplet state consisting of \v e ,i>^), 
-^{\v e ,v e ) + lu^Vfj,)) and \v^v e ), and a singlet state 

-j^{\v e , v e )-\vn, Vfj)). Put another way, the energy eigen- 
states of the system are the usual angular momentum 
states | J, m), where J = 0, 1 and m = — J, — J + 1, . . . , J, 
that carry no "magnetic moment." This is perfectly anal- 
ogous to positronium that consists of two spin-i parti- 
cles with opposite magnetic moments. Neither the sin- 
glet nor the triplet state of positronium carries a mag- 
netic moment. However, we can prepare the system in 
a state with magnetic moment, in our case a state like 
\i/ e ,D e ) = ( 1 1 ■ 0) + 1 , 0) ) . Because the energies of the 
singlet and triplet states are split by the amount /z, we 
will obtain oscillations between the |i/ e ,£/ e ) and \v^,V^) 
states with frequency /i. 

Next, we make the system larger (more extensive) 
without changing its intensive properties, i.e., we keep 
ji = \/2GfN/V fixed with TV the number of neutrino- 
antineutrino pairs and V the volume. If the system is 
prepared in a state consisting of ./V v e -v & pairs, will it 
convert to a state of v^-v^ pairs on a similar time scale 
That is, will the magnetic moment of a "super- 
positronium" consisting of N electrons and N positrons 
reverse on the same time scale as for ordinary positron- 
ium? Formally, this amounts to solving for the expecta- 
tion value (S z — S z ) as a function of time, and, using the 
results of Refs. [3ll.l32l|. it can be shown that the conver- 
sion time scale is of order v^/V/i -1 , not (jT 1 , i.e., much 
longer for a macroscopic ensemble. 

Therefore, quantum effects enter on the usual level of 
1 / y/N fluctuations and do not cause any novel effects on 
a macroscopic scale. Put another way, the equilibration 
in flavour space of a large v e -v e ensemble requires a time 
scale corresponding to ordinary pair processes v e v e «-> 
v^v^ that are of second order in Gf. To first order in Gf, 
the flavour equilibration requires vacuum mixing and the 
phenomenon of bipolar oscillations. 



VIII. DISCUSSION AND SUMMARY 

We have studied bipolar neutrino oscillations, i.e., the 
flavour evolution of an ensemble initially consisting of 
equal numbers of v e and v e . We have shown that the 
classical equations of motion can be cast in the form 



4 References l3lll32H consider a neutrino-only system, initially pre- 
pared with N v £ and M v x , where v x is some linear combination 
of Ve and v^. However, because of the exact correspondence be- 
tween this system and our neutrino-antineutrino system (in the 
absence of Ho) discussed earlier, the results of Refs. [3ll l32f| can 
be trivially mapped to our case. In particular, the connection 
between our (S z — S z ) and their Pi(t) can be found in Sec. HID 
ofRef. [H. 



of a pendulum in flavour space. The surprising bipo- 
lar conversion effect observed for the inverted mass hi- 
erarchy corresponds to the pendulum starting in a near 
upright position, with the excursion angle growing expo- 
nentially until the pendulum makes an almost complete 
swing. Conversely, if it starts in a nearly vertical down- 
ward position (i.e., normal mass hierarchy), the system 
behaves as a harmonic oscillator. 

For the inverted case, the time for a complete swing is 
given by the pendulum's oscillation period times a factor 
depending on the logarithm of the initial excursion angle 
which is nearly identical with twice the vacuum mixing 
angle. Therefore, the bipolar conversion is delayed by 
the logarithm of the small vacuum mixing angle. Like- 
wise, we derived an analytic solution for the case when 
ordinary matter is present and showed that it affects the 
bipolar conversion time also only logarithmically. 

If the vacuum oscillation frequencies are different for 
different modes, one cannot represent the ensemble by 
two "block spins." However, in a dense v-v gas our 
model remains unaffected except that the vacuum os- 
cillation frequency is replaced by an average over all 
modes. When the coupling strength between different 
modes varies, as would be realistically expected in a 
non-isotropic ensemble, yet other forms of behaviour ap- 
pear. In particular, the different modes can kinemat- 
ically decohere in flavour space. Indeed, our results 
suggest — and it has been numerically observed in simula- 
tions [23|, [24| — that partial flavour decoherence, instead 
of a simple swapping of flavours, is a possible feature in 
a multi-coupling/multi-angle system. 

However, the same simulations also suggest that the 
flavour evolution of neutrinos streaming off a supernova 
core is qualitatively approximated by a single-angle treat- 
ment, and that collective behaviour appears to be the 
more generic outcome. A partial explanation for this, at 
least in the inverted hierarchy case, is provided by our ob- 
servation that the single-angle system with a decreasing 
neutrino density and a non-zero neutrino-antineutrino 
asymmetry never becomes properly bipolar. Rather, it 
evolves such that the potential and kinetic energy of 
the pendulum remain equipartitioned, signifying that 
the system remains at the edge of the bipolar condition 
throughout its evolution and is hence better immuncd to 
decoherence. 

The only apparent case of practical interest for this 
discussion is flavour conversion of neutrinos streaming 
off a supernova core where collective flavour transfor- 
mations play an important role. Close to the neutrino 
sphere, the oscillations are synchronised up to a few 
tens of kilometres, then they enter the bipolar regime, 
and finally, beyond 100-200 km, ordinary oscillations oc- 
cur [22|, 123, [2J] • Ordinary matter effects modify the 
oscillations in the usual way both in the synchronised 
regime and far away where collective effects are irrele- 
vant, whereas in the intermediate regime of bipolar oscil- 
lations ordinary matter has no significant impact. This 
counter-intuitive situation was conjectured and numeri- 
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cally observed in Refs. [12, HH, Hit- In our model of a 
flavour pendulum the impact of ordinary matter can be 
calculated analytically. 

Note that while the bipolar behaviour extends over a 
large range in radius outside the neutrino sphere, we have 
explicitly assumed in our treatment that the system con- 
sists of only two flavours because only one mass split- 
ting is of importance. However, it could well be that 
the solar mass difference Am s 2 ok ~ 8x 10 _5 eV 2 can- 
not be ignored in the whole region. If this is the case a 
full three-flavour description must be employed, and new 
phenomena might arise. 

In any case, collective neutrino oscillations, unsup- 
pressed by ordinary matter, in the region a few tens of 
kilometres above the neutrino sphere will likely change 
the picture of supernova flavour oscillations and observ- 
able consequences in various ways. Taking the atmo- 
spheric mass hierarchy to be inverted, the "single-angle 
approximation" predicts a swapping of the D e and 
as well as of the v e and fluxes with a p ossible im- 
pact on r-process nucleosynthesis [26l l27l . l28l [36j , energy 
transfer to the stalling shock wave [371 ]. and the possi- 
bility to observe shock-wave propagation effects in neu- 
trinos [H, [H, EO, SH, S3, [H, [tg, . Nothing new 
happens in the case of a normal mass hierarchy, so that 
one still expects observable effects such as Earth mat- 
ter effects in the neutrino sig nal from the next galactic 
supernova [It], IH, [49|, [5(1 HJ- However, this conclusion 
assumes the validity of the single- angle treatment. Par- 
tial or complete kinematical decoherence, caused by the 
multi-modal nature of the system, will affect the flavour 
composition of the neutrinos passing the bipolar region 
even in the normal hierarchy. 

The one case that probably remains unaffected is the 
prompt deleptonisation burst where initially the v e flux 
is strongly enhanced relative to v u . D, L , v r and v Tl while 
the D e flux is strongly suppressed [52|]. In this case, the 
bipolar condition is not fulfilled and one expects "ordi- 
nary" synchronised oscillations. 
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APPENDIX A: GENERAL EQUATIONS OF 
MOTION 

1. Multiflavour system 

We summarise here the general equations of motion for 
the flavour evolution of an ensemble of mixed neutrinos. 
Our main purpose is to show the meaning of the different 
terms in the general context and their relative signs and 
to establish our conventions. 

A statistical ensemble of unmixed neutrinos is charac- 
terised by the occupation numbers f p = (et p a p ) for each 
momentum mode p, where a p and a p are the relevant cre- 
ation and annihilation operators and (. . .) is the expecta- 
tion value. A corresponding expression can be defined for 
the antineutrinos, f p — (a p a p ), where overbarred quanti- 
ties always refer to antiparticles. In a multiflavour system 
of mixed neutrinos, the occupation numbers are gener- 
alised to density matrices in flavour space [53|, [Hi], 

(p P )ij = {aja^p, 

{Pvh = (a]ai) p . (Al) 

The reversed order of the flavour indices i and j in the 
r.h.s. for antineutrinos is crucial to ensure that p p and p p 
behave consistently under a flavour transformation. The 
seemingly intuitive equal order of flavour indices that is 
frequently used in the literature 

@, SEE!, MM, M, 

EE HI, H2, HH causes havoc in that p q — p* instead of 
Pq — pq appears in Eq. (|A2|) . Therefore, the equations of 
motion then involve p p , p p , /?* and p* and thus lose much 
of their simplicity even if they are, of course, equivalent. 

Flavour oscillations of an ensemble of neutrinos and 
antineutrinos are described by 



dtPp = -i 
dtPp = +i 



O d + \/2Gf \L-L 



d 3 q 
(2tt) 3 
d 3 q 



(Pq - Pq) (! - COs6» pq ) , Pp 



fl p -V2G F (l-L + J ^^(p q -Pq)(l-cos(9 pq )^ , Pp 

r 



(A2) 



where [•,•] is a commutator and Gf is the Fermi constant. cillation frequencies is il p = diag(m 2 , tttJ, m|)/2p with 
For ultrarelativistic neutrinos, the matrix of vacuum os- p — |p| when expressed in the mass basis. The ordinary 
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matter effect is encapsulated in the matrix of charged 
lcpton densities, L — diag(n e , n^, n T ) in the weak inter- 
action basis, and in a corresponding matrix L for the 
charged antilepton densities. The factor (1 — cos6> pq ), 
where # pq is the angle between p and q, will not average 
to unity if the neutrino gas is not isotropic. 

These and the more general Boltzmann kinetic equa- 
tions apply only if no correlations build up between the 
different modes [54] • This condition may well be violated 
when neutrino-neutrino interactions dominate [35| . but 
does not seem to be important in practice for ensembles 
of large numbers of neutrinos 

2. Two-flavour system 

Collective oscillation effects have been studied for the 
case of two-flavour oscillations. The measured hierar- 
chy of neutrino mass differences suggests that oscilla- 
tions driven by the "atmospheric" and "solar" mass dif- 
ferences occur at vastly different epochs in the early uni- 
verse and at vastly different distances from a supernova 
core. Genuine three-flavour collective effects have not 
been addressed in the literature. 

The two-flavour system has the great advantage that 
all 2x2 matrices can be expressed in terms of the unit 
matrix and the Pauli matrices with a vector of coeffi- 
cients. Explicitly we write 

n p = i(fi +w p B-<r) , 
L — L = \ (no + n e L • <r) , 

Pp = M/p + Pp-"")' 

Pp = M/p + Pp'' 7 )- (A3) 



The vectors P p and P p are the v and v polarisation vec- 
tors in flavour space. We choose the coordinate system in 
flavour space such that a polarisation vector pointing in 
the positive z-direction signifies pure electron neutrinos 
or antineutrinos, whereas an orientation in the negative 
z-direction corresponds to muon neutrinos. In this con- 
vention the polarisation vectors are not unit vectors in 
a flavour-pure system. For example, in an ensemble of 
pure electron neutrinos, the z-component of P p corre- 
sponds to the electron neutrino occupation number, and 
the total number density of electron neutrinos would be 
ti„ e = J Pp d 3 p/ (2-7r) 3 . Of course, P p = does not mean 
that this mode is empty; it just means that it contains an 
incoherent equal mixture of electron and muon neutrinos. 

In an ordinary medium there are no charged muons or 
tau-lcptons. Therefore, L is a unit vector in the positive 
z-direction and n e is an effective electron density, i.e., the 
density of electrons minus that of positrons. Finally, vac- 
uum oscillations are determined by the mass differences 
and vacuum mixing angle 6, so that 

w p = (ml — m 2 )/2p, 

B = (sin 2(9,0, cos 260 • (A4) 
Of course, we could have oriented B in any other direc- 
tion in the s-y-plane, i.e., it is our choice to set B y = 0. 
For the normal hierarchy where mi < m,2 the oscillation 
frequency is negative. In the main text we prefer to keep 
a positive ui which implies that we have to reverse the 
z-component of B for the normal hierarchy. 

The terms proportional to the unit matrix in Eq. (|A3|) 
disappear from the equation of motion Eq. (|A2|) due to 
its commutator structure, leaving us with the well-known 
spin-precession equations 



d t P P 
d t P P 



Jp B + V2G F 
j p B -V2G F 



n e L 
rt„L 



dh 



(2tt) 3 
rf 3 q 



(P q -P q )(l-COS0 pq ) 

(P q -P q )(l-COS0 pq ) 



X P r 



X P r 



(A5) 



In the main text we use the frequency 

A = V2G F n e (A6) 

as a coefficient for L to quantify the matter effect. An 
ensemble consisting initially of v e and D e corresponds to 
J P p tf 3 p/(2-7r) 3 = (0, Q,n Ve ). Therefore, if we represent 
the entire v e ensemble with a single integrated polarisa- 
tion vector P of unit length, the v-v term must be of the 
form fi(P — P) x P p , with 

fi = V2G F n„ c = V2G F 7i Pe . (A7) 

We use this frequency to denote the strength of the neu- 
trino self coupling. 



The equations of motion of the entire system, assuming 
all neutrinos have the same vacuum oscillation frequency, 
thus become 

d t V = [+cjB + AL + /i (P - P)] x P, 

d t P = [— cjB + AL + /i (P - P)] x P . (A8) 

The B-parts of the Hamiltonian and of the equations of 
motion correspond exactly to those of a particle and its 
antiparticle with a magnetic moment in the presence of a 
B-field. They have opposite magnetic moments and thus 
spin-precess in opposite directions. 

The matter term includes an important sign-change in 
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that the particle and antiparticle have equal energies if 
their spins are aligned, but their magnetic moments are 
anti-aligned. Of course, this sign change reflects that in 
the presence of a medium, particles and antiparticles are 
affected in opposite manners relative to the vacuum term 
so that the usual MSW effect occurs for the normal, but 
not the inverted mass hierarchy. 



3. Supernova Neutrinos 

Bipolar oscillations are primarily important for neu- 
trinos streaming off a supernova core. Therefore, we 
briefly state the typical parameter values expected in this 
context. In numerical simulations of supernova neutrino 
oscillations, it is often assumed that (E Ve ) = 11 MeV, 
(Eu e ) = 16 MeV and {E Vx ) = 25 MeV^for the other 
species [23|, [29| . The "atmospheric" neutrino mass differ- 
ence relevant here is Am| 3 = 1.9-3.0 x 10~ 3 eV 2 . With 



1.9-3.0 x 10" 

(E v ) = 16 MeV, we may thus use 

0.3 km" 1 



(A9) 



as a typical number. In a supernova one studies the neu- 
trino flavour evolution as a function of radius from the 
neutrino sphere, so it is useful to express all distances in 
km and all frequencies in km -1 . 

Moreover, most numerical simulations assume that all 
neutrino species are emitted with the same luminosities, 
with Lq = 10 51 erg s _1 being a typical choice. The 
neutrino-sphere radius is approximately R u — 10 km. 
Therefore, a typical neutrino density at radius r is 
n v = L /({Ej)4nr 2 ) = 1.04 x 10 32 cm' 3 ^ 2 , where 
rio = r/R u = r/10 km. The relevant density for the 
calculation of fi is the difference between the v e and the 
densities. This amounts to reducing n v by a factor 



16/11 — 16/25 f» 0.81, a number that reflects the different 
average energies of the different species. Finally, we need 
to include the typical angular factor 1 — cos 6 pq between 
neutrino trajectories because collinear neutrinos do not 
cause refractive effects for each other. This angular ef- 
fect is approximately taken into account with the factor 
F = |[1 — (1 — Rl/r 2 ) 1 / 2 ] 2 used in previous "one-angle" 
numerical studies and originally worked out in Ref . [27[ . 
Altogether we thus find that 



H = 0.3 x 10 5 km -1 1 - (1 - 



M 2 ' 



(A10) 



is a reasonable value for simple estimates. Near the neu- 
trino sphere, [i is 10 5 times larger than to. 

The neutrinos streaming off a supernova core are not 
isotropic so that the "multi-angle" nature of the problem 
can be important. Still, in a given radial direction, the 
problem can be assumed to have axial symmetry so that 
different neutrino modes can be classified by their angle 
Cj = cos 6i , where is the angle of the neutrino mo- 
menta relative to the radial direction, i.e., c; represents 
all neutrinos streaming in the direction 9i, integrated over 
all azimuthal directions. The coupling strength for two 
different "modes" 9i and 9j , weig hted by the angular fac- 
tors, is then proportional to [231 ] 



1 CiCj 



(All) 



Of course, in an isotropic medium, where Ci and Cj are 
both uniformly distributed between —1 and +1, this term 
averages to 1. However, if we consider neutrinos emitted 
isotropically from a flat surface, we will have a uniform 
distribution in the range < Cj < 1. This provides a 
simple model for a nonisotropic medium. 



[1] L. Wolfenstein, "Neutrino oscillations in matter," Phys. 
Rev. D 17, 2369 (1978). 

[2] S. P. Mikheev and A. Y. Smirnov, "Resonance enhance- 
ment of oscillations in matter and solar neutrino spec- 
troscopy," Yad. Fiz. 42, 1441 (1985) [Sov. J. Nucl. Phys. 
42, 913 (1985)]. 

[3] S. P. Mikheev and A. Y. Smirnov, "Resonant amplifica- 
tion of neutrino oscillations in matter and solar neutrino 
spectroscopy," Nuovo Cim. C 9, 17 (1986). 

[4] G. M. Fuller, R. W. Mayle, J. R. Wilson and D. N. 
Schramm, "Resonant neutrino oscillations and stellar col- 
lapse" Astrophys. J. 322 (1987) 795. 

[5] D. Notzold and G. Raffelt, "Neutrino dispersion at finite 
temperature and density," Nucl. Phys. B 307, 924 (1988). 

[6] T. K. Kuo and J. T. Pantaleone, "Neutrino oscillations in 
matter," Rev. Mod. Phys. 61, 937 (1989). 

[7] J. Pantaleone, "Neutrino oscillations at high densities," 
Phys. Lett. B 287, 128 (1992). 

[8] S. Samuel, "Neutrino oscillations in dense neutrino gases," 
Phys. Rev. D 48, 1462 (1993). 



[9] V. A. Kostelecky, J. Pantaleone and S. Samuel, "Neutrino 
oscillation in the early universe," Phys. Lett. B 315, 46 
(1993). 

[10] V. A. Kostelecky and S. Samuel, "Neutrino oscillations 
in the early universe with an inverted neutrino mass hier- 
archy," Phys. Lett. B 318, 127 (1993). 

[11] V. A. Kostelecky and S. Samuel, "Nonlinear neutrino os- 
cillations in the expanding universe," Phys. Rev. D 49, 
1740 (1994). 

[12] V. A. Kostelecky and S. Samuel, "Self-maintained coher- 
ent oscillations in dense neutrino gases," Phys. Rev. D 52, 
621 (1995) [hep-ph/9506262] , 

[13] V. A. Kostelecky and S. Samuel, "Neutrino oscillations 
in the early universe with nonequilibrium neutrino distri- 
butions," Phys. Rev. D 52, 3184 (1995) [hep-ph/9507427] , 

[14] S. Samuel, "Bimodal coherence in dense selfinteract- 
ing neutrino gas es," Phys. Rev. D 53, 5382 (1996) 
[hep-ph/9604341| . 

[15] V. A. Kostelecky and S. Samuel, "Nonequilibrium neu- 
trino oscillations in the early universe with an inverted 



21 



neutrino mass hierarchy," Phys. Lett. B 385, 159 (1996) 
[hep-ph/9610399] . 
[16] J. Pantaleone, "Stability of incoherence in an isotropic 
gas of oscillating neutrinos," Phys. Rev. D 58, 073002 
(1998). 

[17] S. Pastor, G. G. Raffelt and D. V. Semikoz, 
"Physics of synchronized neutrino oscillations caused 
by self-interactions," Phys. Rev. D 65. 053011 (2002) 
[hep-ph/0109"035] . 

[18] C. Lunardini and A. Y. Smirnov, "High-energy neutrino 
conversion and the lepton asymmetry in the universe," 
Phys. Rev. D 64, 073006 (2001) [hep-ph/0012056] . 

[19] A. D. Dolgov, S. H. Hansen, S. Pastor, S. T. Petcov, 
G. G. Raffelt and D. V. Semikoz, "Cosmological bounds 
on neutrino degeneracy improved by flavor oscillations," 
Nucl. Phys. B 632, 363 (2002) [hep-ph/0201287| . 

[20] Y. Y. Y. Wong, "Analytical treatment of neutrino 
asymmetry equilibration from flavour oscillations in 
the early univers e," Phys. Rev. D 66, 025015 (2002) 
[hep-ph/ 0203180] . 

[21] K. N. Abazajian, J. F. Beacom and N. F. Bell, "Stringent 
constraints on cosmological neutrino antineutrino asym- 
metries from synchronized flavor transformation," Phys. 
Rev. D 66, 013008 (2002) [astro-ph/0203442] . 

[22] H. Duan, G. M. Fuller and Y. Z. Qian, "Col- 
lective neutrino flavor transformation in supernovae," 
astro-ph/0511275| 

[23 



H. Duan, G. M. Fuller, J. Carlson and Y. Z. Qian, "Sim- 
ulation of coherent non-linear neutrino flavor transforma- 
tion in the supernova environment. I: Correlated neutrino 
trajectories," astro-ph/0606616 
[24] H. Duan, G. M. Fuller, J. Carlson and Y. Z. Qian, "Co- 
herent development of neutrino flavor in the supernova 
environment," astro- ph/0608050 
[25] M. T. Keil, G. G. Raffelt and H. T. Janka, "Monte Carlo 
study of supernova neutrino spectra formation," Astro- 
phys. J. 590, 971 (2003) [astro-ph/0208035] . 
[26] J. T. Pantaleone, "Neutrino flavor evolution near a 
supernova's core," Phys. Lett. B 342, 250 (1995) 
[astro-ph/9405008] . 
[27] Y. Z. Qian and G. M. Fuller, "Neutrino-neutrino scat- 
tering and matter enhanced neutrino flavor transforma- 
tion in Supernova e," Phys. Rev. D 51, 1479 (1995) 
[astro-ph/9406073] . 
[28] G. Sigl, "Neutrino 
nova nucleosynthesis," 
[astro-ph/9410094] 
[29] S. Pastor and G. Raffelt. "Flavor oscillations in the 
supernova hot bubble region: Nonlinear effects of neu- 
trino background," Phys. Rev. Lett. 89, 191101 (2002) 
[astro-ph/0207281] . 
[30] A. B. Balantekin and H. Yiiksel, "Neutrino mixing and 
nucleosynthesis in core-collapse supernovae," New J. Phys. 
7, 51 (2005) [astro-ph/0411159] . 
[31] A. Friedland and C. Lunardini, "Do many-particle neu- 
trino interactions cause a novel coherent effect?," JHEP 
0310, 043 (2003) [hep-ph/0307140] . 
[32] A. Friedland, B. H. J. McKellar and I. Okuniewicz, 
"Construction and analysis of a simplified many-body 
neutrino model," Phys. Rev. D 73, 093002 (2006) 
[hep-ph/0602"016] . 
[33] A. BT Balantekin and Y. Pehlivan, "Neutrino neu- 
trino interactions and flavor mixing in dense matter," 
|astro-ph/0607527| 



mixing constraints 
Phys. Rev. D 51, 



and 
4035 



super- 
(1995) 



[34] A. Friedland and C. Lunardini, "Neutrino flavor con- 
version in a neutrino background: Single- versus multi- 
particle descripti on," Phys. Rev. D 68, 013007 (2003) 
[hep-ph/0 304055 . 
[35] N. F. Bell, A. A. Rawlinson and R. F. Sawyer, "Speed-up 
through entanglement: Many-body effects in neutrino pro- 
cesses," Phys. Lett. B 573, 86 (2003) [hep-ph/0304082] . 
[36] Y. Z. Qian, G. M. Fuller, G. J. Mathews, R. Mayle, 
J. R. Wilson and S. E. Woosley, "A connection between 
flavor mixing of cosmologically significant neutrinos and 
heavy element nucleosynthesis in supernovae," Phys. Rev. 
Lett. 71, 1965 (1993). 
[37] G. M. Fuller, R. Mayle, B. S. Meyer and J. R. Wilson, 
"Can a closure mass neutrino help solve the supernova 
shock reheating problem?", Astrophys. J. 389, 517 (1992). 
[38] R. C. Schirato and G. M. Fuller, "Connection between 
supernova shocks, flavor transformation, and the neutrino 
signal," |astro-ph/0205390 



E. Dalhed and J. R. Wil- 
and neutrino oscillation 
Phys. 20, 189 (2003) 



[39] K. Takahashi, K. Sato, H 
son, "Shock propagation 
in supernova," Astropart 
as tro-ph/0212195] . 

[40] C. Lunardini and A. Y. Smirnov, "Probing the neu- 
trino mass hierarchy and the 13-mixing with supernovae," 
JCAP 0306, 009 (2003) [hep-ph/0302033] . 

[41] G. L. Fogli, E. Lisi, D. Montanino and A. Mirizzi, "Anal- 
ysis of energy- and time-dependence of supernova shock 
effects on neutrino crossing probabilities," Phys. Rev. D 
68, 033005 (2003) [hep-ph/0304"056] 

[42] R. Tomas, M. Kachelriefi, G. Raffelt, A. Dighe, H. T. 
Janka and L. Scheck, "Neutrino signatures of supernova 
shock and reverse shock propagation," JCAP 0409, 015 
(2004) [astro-ph/0407132] . 

[43] B. Dasgupta and A. Dighe, "Phase effects in neu- 
trino conversions during a supernova shock wave," 
|hep-ph/0510219| 

[44] G. L. Fogli, E. Lisi, A. Mirizzi and D. Montanino, "Prob- 
ing supernova shock waves and neutrino flavor transitions 
in next-generation water-Cherenkov detectors," JCAP 
0504, 002 (2005) [hep-ph/0412"046] , 

[45] S. Choubey, N. P. Harries and G. G. Ross, "Probing neu- 
trino oscillations from supernovae shock waves via the Ice- 
Cube detector," |hep-ph /0605255 

[46] G. L. Fogli, E. Lisi, A. Mirizzi and D. Montanino, 
"Damping of supernova neutrino transitions in stochas- 
tic shock-wave de nsity profiles," JCAP 0606, 012 (2006) 
[hep-ph/0603033l . 

[47] A. S. Dighe and A. Y. Smirnov, "Identifying the neutrino 
mass spectrum from the neutrino burst from a supernova," 
Phys. Rev. D 62, 033007 (2000) [hep-ph/9907423) . 

[48] C. Lunardini and A. Y. Smirnov, "Supernova neutrinos: 
Earth matter effects and neutrino mass spectrum," Nucl. 
Phys. B 616, 307 (2001) [hep-ph/0106"l49] . 

[49] A. S. Dighe, M. T. Keil and G. G. Raffelt, "Detecting the 
neutrino mass hierarchy with a supernova at IceCube," 
JCAP 0306, 005 (2003) [hep-ph/0303210] . 

[50] S. H. Chiu and T. K. Kuo, "Probing neutrino mass hier- 
archies and 0i3 with supernova neutrinos," Phys. Rev. D 
73, 033007 (2006); Erratum ibid. D 73, 059901 (2006). 

[51] A. Mirizzi, G. G. Raffelt and P. D. Serpico, "Earth mat- 
ter effects in supernova neutrinos: Optimal detector loca- 
tions," JCAP 0605, 012 (2006) [astro-ph/0604300| . 

[52] M. Kachelriefi, R. Tomas, R. Buras, H. T. Janka, 
A. Marek and M. Rampp, "Exploiting the neutronization 



22 



burst of a galactic supernova," Phys. Rev. D 71, 063003 

(2005) [astro-ph/0412082| . 
[53] A. D. Dolgov, "Neutrinos in the early universe," Yad. 

Fiz. 33, 1309 (1981) [Sov. J. Nucl. Phys. 33, 700 (1981)]. 
[54] G. Sigl and G. Raffelt, "General kinetic description of rel- 



ativistic mixed neutrinos," Nucl. Phys. B 406, 423 (1993). 
[55] B. H. J. McKellar and M. J. Thomson, "Oscillating dou- 
blet neutrinos in the early universe," Phys. Rev. D 49, 
2710 (1994). 



APPENDIX B: ERRATUM 
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In Eq. (46) of "Self-induced conversion in dense neu- 
trino gases: Pendulum in flavour space" [l[ we give 



(i 



transition 



4(1 + a) 



(Bl) 



for the transition between the synchronized and bipolar 
behavior of the flavor pendulum. Here, u) is the vacuum 
oscillation frequency, /i quantifies the neutrino-neutrino 
interaction strength, and < a < 1 where a = rip/n^. 
The neutrino gas is "dense" when w/n <C 1. Equa- 
tion (|BTj) was derived under this assumption and as such 
is self-consistent only if a is not too small. 

Subsequently it was noted Q that this approximation 
is not necessary. In the relevant case of a vanishing mix- 
ing angle, we have a = 1 — a in Eq. (45) of Ref. [l| and 
Q = (1+a— £) for the inverted hierarchy so that Eq. (45) 
is equivalent to 



(l-a) 2 =4£(l + «-£)• 



It is solved by 



(i - VS) 2 



(B2) 



(B3) 



in agreement with Eq. (76) of Ref. Q. 

To compare Eq. (|B1[) with its alternative form 
Eq. (|B3[) , we use a = 1 — a and expand, 



c = a l_ f l + |o-+icr 2 4 
"~ 8 X \l + |a+^<r 2 



from Eq. (|ET 



from Eq. (JB3 



(B4) 



Both results are surprisingly similar. For the example 
a = 0.8 used in Refs. [H, 0, the approximate result from 
Eq. ([BTj) is 0.3% smaller than the exact one of Eq. (lB3l . 



As noted in Ref. 0, Fig. 8 and the text around Eq. (49) 
of our paper are not consistent with the definitions of the 
energies in Eqs. (47) and (48). In Fig. 8 and its discus- 
sion, we have used E pot = cu B • (P + P) for the potential 
energy, i.e., the first part of Eq. (34), and the second part 
of Eq. (34) for the kinetic energy. In this way the poten- 
tial energy does not depend on /i, and the kinetic energy 
vanishes for /i — > 0. On the other hand, the flavor pendu- 
lum is described by Q = P + P — (oj//i)B. Its potential 
energy is wB • Q. When lo/^l is not small, there is an 
important difference between these descriptions. For a 
system with fixed \x one can add arbitrary functions of /x 
to the potential or to the kinetic energy without affecting 
the dynamics. On the other hand, when /i and thus the 
Hamiltonian depend on time, one has to be more careful. 
The main point of Fig. 8, the discovery of the approxi- 
mate equipartition of energies as the system evolves, re- 
mains unaffected and is analytically elaborated in Ref. 
in terms of adiabatic invariants. 

As noted in Ref. [|[ , another clarification concerns the 
normalization of the energies. Equation (34) and subse- 
quently all expressions for E pot and i?kin should include 
a factor 1/2 if the energy is to be interpreted as an en- 
ergy per neutrino. The quantum Hamiltonian Eq. (67) is 
normalized without ambiguity. The normalization of the 
classical energy is inconsequential for our results. 
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